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Abstract 

Q ■ The problem of structure formation in relativistic dissipative fluids was analyzed in a pre- 

^ ■ vious work within Eckart's framework, in which the heat flux is coupled to the hydrodynamic 

^ acceleration, additional to the usual temperature gradient term. It was shown that in such 

case, the pathological behavior of fluctuations leads to the disapperance of the gravitational 
instability responsible for structure formation In the present work the problem is re- 
visited now using a constitutive equation derived from relativistic kinetic theory. The new 
relation, in which the heat flux is not coupled to the hydrodynamic acceleration, leads to 
a consistent first order in the gradients formalism. In this case the gravitational instability 
remains, and only relativistic corrections to the Jeans wave number are obtained. In the 
calculation here shown the non-relativistc limit is recovered, opposite to what happens in 
Eckart's case 
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I. INTRODUCTION 



Relativistic hydrodynamics has been a subject of interest since the 1940 's. However, 
many open questions remain unanswered, in particular the nature of heat conduction in hot 
gases has recently been under research and debate. The main issue in such a subject concerns 
the particular structure of the constitutive equation for dissipative energy flow in the special 
relativistic case. On the one hand, C. Eckart in 1940, proposed a relation which coupled 
heat to the hydrodynamic acceleration [3]. This relation was obtained phenomenologically 
as a sufficient condition for the entropy production to remain positive. Later in 1985, 
Hiscock and Lindblom working relativistic dissipative fluids, found generic instabilities in 
!5ckart's formalism which in part lead to the developement and use of extended theories 
A subsequent study showed that these instabilities are due to the presence of the 
hydrodynamic acceleration in the heat flux constitutive equation j4|. 

On the other hand several authors have found that kinetic theory predicts a coupling 
of heat with gradients of state variables instead of the acceleration p-lT]. At the present 
time, both type of constitutive equations have been used in order to study the Rayleigh- 
Brillouin scattering spectrum, finding that by introducing Eckart's equation the results 
were inconsistent with the observations, while the use of the equations in terms of gradients 
eliminates the inconsistencies Q|. 

In recent work, two of us addressed the question of structure formation in relativistic 
dilute gases in the presence of dissipation introducing Eckart's equation for the heat fiux 
{l|. The results obtained were similar to those presented by Hiscock & Lindblom, that is 
the apperance of an instability of early onset triggered by the fiuid itself. This instability, 
additional to not being observed, rules out any possibility of structure formation in such 
gases since the growth of density fiuctuations due to the gravitational field are overriden. 

In the present work, we instead use the consitutive equation obtained in relativistic kinetic 
theory in the study of structure formation and find that generic instabilities are eliminated in 
such a way so that gravity becomes dominant. We finally establish the relativistic corrections 
to the critical wave number. 

We have structured this paper as follows. In Sect. II we establish the system of relativistic 
transport equations for a gas including gravitational effects within the Newtonian limit and 
linearize around equilibrium values for the state variables in order to obtain a set of linear 
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equations for the fluctuations. In Sect. Ill we calculate the dispersion relation for the system 
in the Fourier-Laplace space and infer from its solutions the behavior of density fluctuations. 
Section IV is devoted to the explicit calculation of the relativistic correction to the Jeans wave 
number in a dissipative medium and compare the result with its non-relativistic counterpart. 
Conclusions and flnal remarks are included in Sect. V. 



II. LINEARIZED TRANSPORT EQUATIONS 

The behavior of density fluctuations in a fluid is analyzed by studying the response of 
the corresponding linearized system of transport equations when the state variables are 
perturbed from their equilibrum values j^. That is, if X is a state variable one considers 

X = Xo + 5X (1) 

where Xq is the equilibrium value and 5X the fluctuation. The gas here considered is a 
dilute cloud of particles with a temperature high enough, so that the relativistic effects in 
the molecular motion is relevant. We consider the temperature T, number density n and 
hydrodynamic four- velocity as the independent, local equilibrium, state variables. The 
system of equations where the hypothesis in Eq. ([1]) is to be introduced is composed of two 
tensor conservation equations namely, 

K = (2) 

and 

= (3) 

where 

N" = nu^ (4) 

is the particle four-flow and 

TIF 11 

= —u^u, + pK + < + -oQ^u, + -u^'q, (5) 

is the stress-energy tensor. Here e is the energy density per particle, c the speed of light and 
p the hydrostatic pressure. The dissipative fluxes are the Navier tensor t^^'^ and the heat 
flux q^. We have also used the standard expression for the spatial projector 
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in a + + +— metric and the orthogonality relations 



u'^qy = u^TCf,,, = 
are also assumed. Thus, the transport equations for the relativistic fluid may be written as 

fi + ne = (6) 



TIE V \ ( ^E J) \ 

^ + ^ j + + f,e) Uy + P,,K + 

+ \ + (l^u,,,^ + + W'qu,^) = (7) 



nCnf + p9 + u-f^TTi: + + ^u'q, = (8) 

where the relation e = CnT has been used in order to turn the internal energy balance 
equation in the heat equation. Also, we have deflned 6 = u^^ and a dot denotes a proper 



equations, constitutive relations for the 
I 

1| the consitutive equations proposed by 



time derivative. In order to complete the set o 
dissipative terms have to be provided. In Ref. 
Eckart, namely 

< = -2r^KKr^ - (66^, (9) 

and 

= -^h; (t'^ + ^^u^^ (10) 

where introduced. In the present work we instead use the constitutive relation obtained from 
kinetic theory by solving Botlzmann's equation by the Chapman- Enskog approximation, that 
is 

q" = -h';(LTT— + L^T—) (11) 



a. 



for the heat flux, while the equation for the momentum flux identical 

In order to account for self-gravitational effects prior to collapse, we consider a Newtonian 
approximation of general relativity given by the line element 

ds' = dr^ + r'de^ + sin' ed(j)' - (^1 - {cdtf (12) 

where ip is the ordinary graviational potential. Since, to flrst order in fluctuations, the con- 
tinuity and heat equations, Eqs. and (IHl), are only coupled with the velocity through its 



divergence 59, it is convenient to calculate the divergence and curl of the momentum balance 
equation Eq. ([7]). It can be shown that by performing such calculation, the equation for 
the curl of the velocity decouples from the system and has decaying solutions for the corre- 
sponding fluctuations On the other hand, the equation for 66 now has an inhomogeneous 
term that can expressed in terms of density fluctuations by using the Poisson equation. The 
resulting linearized system of equations reads 

^(<5n)+«o^^ = (13) 

Po|- ((5^)+fcToV' {6n)+nQkV^ {6T)-AV^ {Se)-LTTV^ (6t)-^V^ {6n) = -AwGrnpo {5n) 

(14) 

ano^ (5T) + nokTo (66) - LttV' {6T) - L^t^^ {5n) = (15) 

where, for convenience we have deflned po = {nQEQ +po) /c^ and considered a fluid at rest, 
or an observer in a comoving frame, such that the equilibrium value of the hydrodynamic 
velocity vanishes. Equations fll3tll5p constitute the system of linearized equations for the 
relativistic fluid to flrst order in the gradients considering a weak (Newtonian) gravitational 
fleld. In the next section, the set will be analyzed in Laplace-Fourier space in order to 
address the behavior of density fluctuations. 



III. DISPERSION RELATION 

In the previous section, the system of transport equations to flrst order, both in gradients 
and in fluctuations has been established. In order to analyze the behavior of_density fluc- 
tuations, we proceed by following the standard method in fluctuation theory [9]. The set is 
transformed to Laplace-Fourier space and the dispersion relation is obtained by setting the 
corresponding determinant equal to zero. In this case, such a relation is given by 

/ . ^ 



det 



no 



kToq^ + ^sq^ + AnGmpo Aq^ + pos -nokq^ 



-sq- 



V 



nokTo 







(16) 



\knoS + LttQ^ J 



or, as a cubic equation 



+ a2 (g) + ai (g) s + oq (g) = 



(17) 
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where the coefficients are given by 



and 



02 g = — + —p^ Q (18) 

PoC^C„ 3po 

ao (g) = g q (20) 

It is worthwhile pointing out at this stage that by considering A = L^t = L^t = for the 
non-dissipative case, one is led directly to the Jeans wave number 

d = (21) 

where 

is the speed of sound in the medium, such that the only relativistic corrections to the 
Jeans instability criterion neglecting dissipation arise from the relativistic values of po and 
Cn- Indeed, as is verifiable from the system of equations, the relativistic effects are mostly 
coupled with the dissipative ones. Since in the non-relativistic limit po ~ mriQ we recover 
the usual Jeans wave number, where Ct = W^- 

Turning back to the complete dispersion relation, one can extract information about the 
roots of such equation by using the results obtained in the case without gravitational field. 
We recall that in such calculation, both relativistic and non-relativistic dispersion relations 
have one real solution, which gives rise to a central peak called the Rayleigh peak, and two 
imaginary roots which lead to damped acustic modes corresponding to the Brillouin peaks. 
For a more detailed analysis of these solutions and the corresponding scattering spectrum 
the reader may consult Refs. jsl, jol and 



10|. Based on those results, we here assume that 



the Rayleigh peak is not significantly modified by the presence of the gravitational field such 
that one of the roots of Eq. f lTTl) still has the form s = jq'^ and find that this condition is 
satisfied by the coefficient 

^ ^ Ltt ^23) 

n^Cn 

Under this assumption, we can factorize Eq. (fTTl) as follows 

(s + 7g2) (s^ + ps + i^) = (24) 



and, keeping only up to second order in q we obtain 

PO V C2 C^Cn 



= ^ M - -5 37^ ? (25) 



and 



u = -ATcGmno + ^ (l + ^]q' (26) 



The factorization in Eq. fl24p corresponds to the original dispersion relation given in Eq. 
( IT7|) up to order which is consistent with the approximation here considered. Having two 
roots of the equation, arising from a second order polynomial yields the same physics as in 
the Jeans instabiliy with and without dissipation, since one can find limiting values for q for 
which either oscillating or growing modes exist. 

The limiting value for growing modes in density fluctations, which physically leads to 
a gravitational collapse, is extracted directly from the discriminant of the second order 
polynomial in Eq. (fT7|) . that is 



'a - ^ - ' q' - kT,n,~p, U + + Ar.Gmn,~pl = (27) 



which has a root given by 



47rG' ^ / 2noLnT ^TqLtt 



(28) 



As mentioned before, this wave number corresponds to the limit between oscillating and 
exponentially gowing modes. At this point it is convenient to point out that in previous 
work the dispersion relation obtained was a fourth order polynomial with a dominating 
positive root which wiped out any posibility of structure formation in the system. On the 
other hand, the use of Eq. (fTT]) leads to the same behavior as in the non-relativistic case, 
with small relativistic corrections as will be shown in the next section. 



IV. CORRECTION TO THE JEANS WAVE NUMBER AND NON- 
RELATIVISTIC LIMIT 



In order to assess the magnitude of the corrections to the Jeans wave number implied by 
Eq. (128|) we expand the terms in the square root, so that 



nGA ( ^ _ 2noLnT ^TqLtt 



CtPo 



3c2 



(29) 
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where 



~2 47rGpo /o^^ 



is the Jeans wave number corrected to the relativistic values of po and C]. Notice that the 
corrections beyond ~ qj arise solely from the coupling of the gravitational effects with the 
dissipative ones. Moreover, the relativistic corrections depend of the heat flux but are also 
coupled with the viscous effect. This is an interesting result which gives further insight in 
the role of viscous dissipation in relativistic systems and reassures the need of more research 
in this direction. 

The non-relativistic limit is readily obtained when c — )■ oo since — > qj and we obtain 



TT 



"^^'H^-^j '''' 

which is precisely the result in Ref. ll| with = qj included in Eq. (13 ip . for the case of 
a non-dissipative medium. Thus, the introduction of a constitutive equation for the heat 
flux in terms of gradients of scalar state variables - temperature, density, pressure - assures 
that the gravitational collapse is possible for relativistic systems following roughly the same 
physics as in the non-relativistic scenario, that is with the existence of a critical wavenumber 
indicating a range of values for which a structure will form within the gas. 



V. CONCLUSIONS AND FINAL REMARKS 

The issue of which is the correct formalism for describing dissipative processes in rela- 
tivistic hydrodynamics has recently been a subject of intense research. It was already shown 
that the causality and stability issues allegedly present in first order theories are not a con- 
cern when the constitutive equations arising from the microscopic formalism are considered 

13|. The causality of the system is guaranteed as long as one considers fluctuations in 

I, 

the hydrodynamic velocity even in a comoving frame This is a completely reasonable 
assumption since either a fluid at rest or an observer in a comoving frame shall be stated as 
a vanishing equilibrium velocity. Fluctuations, which are physical evidence of the statistical 
picture underlying the hydrodynamic formalism, are spontaneous and always present. Thus, 
there is really no reason to consider 66 = 0. On the other hand, the stability concern is 
actually more a problem about the relaxation of such fluctuations instead of a hydrodynamic 
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instability per-se. In any case, it has been shown that the pathological behavior of the sys- 
tem is wiped out when gradients are included in the momentum balance through the heat 
lux constitutive equation instead of the acceleration which appears in Eckart's formalism 
^. These facts lead to the firm statement that extended theories are not strictly neccessary 
and first order theories are physically sound as long as the results of kinetic theory are used 
to close the system of transport equations. However, as shown in Ref. jl2|, introducing 
relaxation parameters also solves the problem by forcing into the equations a decay that 
overrides the instability. 

The question of which formalism is more adeccuate is still a subject of debate. In the 
present paper the calculation presented in Ref. |l| by two of us is carried out once again, but 
now using the closure relation obtained from kinetic theory in order to argue in favor of first 
order theories. It is a fact that structures form in systems when the Jeans criterion is met. 
It is desirable to have also this kind of behavior in the relativistic formalism with corrections 
which vanish for mild or low temperatures. In Ref. jl| it was clearly shown that if Eckart's 
closure is introduced in the transport equations, the gravitational instability is not present 



in the system, even in the non-relativstic limit. On the other hand, wha'; 



here is that by introducing the constitutive equation obtained in Refs. 



las been shown 



471 , the physics in 



the gravitational collapse in a relativistic system matches the phenomenon established by 
Jeans in the non-dissipative case jl5| as well as the one described by Sandoval- Villalbazo and 
Garcfa-Colfn for dissipative systems The mechanism for the onset of the gravitational 
instability is the same and only slight corrections to the Jeans wave number are obtained. 
It is of the authors' opinion that even if these corrections are not astrophysically significant, 
the fact that the theory predicts a collapse under particular conditions and that it yields 
the correct non-relativistic limit is enough evidence to rule it as more adequate to describe 
the physics of relativistic gases than other proposals. 
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